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Abstract 



The processes e + e~ — > 777, 7r°7 — > 37 have been studied in the cm. energy range 
600-1380 MeV with the CMD-2 detector. The following branching ratios have been 
determined: 



B(p° -► 777) = (3.21 ± 1.39 ± 0.20) • 10 -4 , 
B{uj -». r? 7 ) = (4.44±?;|| ± 0.28) • 10~ 4 , 

-► 777) = (1.373 ± 0.014 ± 0.085) • 10~ 2 , 
B(p° 7r°7) = (6.21+{;^ ± 0.39) • 10" 4 , 
B(w -► vr°7) = (9.06 ± 0.20 ± 0.57) • 10~ 2 , 
B(<j> -► tt°7) = (1.258 ± 0.037 ± 0.077) • 10~ 3 . 



1 Introduction 



The magnetic dipole transitions of the light vector mesons (p, u and 0) to 
the 7r°7 and 777 final states have traditionally provided a convenient labora- 
tory for various tests of theoretical concepts, particularly the nonrelativistic 
quark model and Vector Dominance Model (VDM) [1,2]. There are ongoing 
discussions about mechanisms of SU(3) breaking, possible admixture of glue 
in mesons and the role of anomalies in radiative decays [3,4,5,6,7,8]. Precise 
measurements of the cross sections of e + e~ annihilation into the 7r°7 and 777 
final states in the broad c.m.energy range are necessary for the problem of the 
muon anomaly [9] . Radiative decays to 7r°7 and 777 can also provide important 
information on the properties of the p, u and excitations as well as on the 
existence of light hybrids between 1000 and 2000 MeV [10,11]. 



Despite previous experimental efforts (cf. the detailed bibliography in [12]), 
of these decays only 00 — > 7r°7 and — > 777 are rather well studied. A three- 
photon final state is convenient for the investigation of the 7r°7 and 777 final 
states since both 7r° and rj readily decay into two photons. Measurements of 
the branching ratios for corresponding decays of the p, u and <fi using the 
two-photon decay mode have been performed at ND [13,14] and SND [15,16], 
however, none of them covered the whole off-resonance energy range. 



In this work we report on the measurement of the cross section of the processes 
e + e~ — > 7r°7 and e + e~ — > 777 in the three-photon final state in the c.m.energy 
range 600-1380 MeV using the data from the CMD-2 detector at the VEPP- 
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2M e + e" collider. 



2 Experiment 

The general purpose detector CMD-2 has been described in detail elsewhere [17] 
Its tracking system consists of a cylindrical drift chamber (DC) and double- 
layer multiwire proportional Z-chamber, both also used for the trigger, and 
both inside a thin (0.38 X ) superconducting solenoid with a field of 1 T. 
The barrel Csl calorimeter (BC) with a thickness of 8.1 X placed outside the 
solenoid has energy resolution for photons of about 9% in the energy range 
from 100 to 700 MeV. The angular resolution is of the order of 0.02 radians. 
The end-cap BGO calorimeter with a thickness of 13.4 X placed inside the 
solenoid has energy and angular resolution varying from 9% to 4% and from 
0.03 to 0.02 radians, respectively, for the photon energy in the range 100 to 
700 MeV. The barrel and end-cap calorimeter systems cover a solid angle of 
0.92 x 47r radians. 

This analysis is based on a data sample corresponding to integrated luminos- 
ity of 21 pb" 1 collected in 1997-1998 in the energy range 600-1380 MeV. The 
step of the cm. energy scan varied from 0.5 MeV near the oo and <fi peaks to 
10 MeV far from the resonances. The beam energy spread is about 4 x 10~ 4 
of the total energy. The luminosity is measured using events of Bhabha scat- 
tering at large angles [18]. 

A GEANT3 based Monte Carlo simulation (MC) package is used to model the 
detector response and determine the efficiency [19]. Because of the beam in- 
duced background additional ("fake") clusters can appear in the calorimeter. 
To take this effect into account in MC we determine a corresponding probabil- 
ity as well as photon energy and angular spectra directly from the data using 
the process e + e~ 7r + 7r~7r°, and then include generation of such photons 
in the detector response during simulation. 



3 Data analysis 

At the initial stage, events are selected which have no tracks in the DC, three 
or four photons in the Csl calorimeter, the total energy deposition 0.8 < 
E tot /E cm < 1.1, the total momentum P tot /E cm < 0.15 and the minimum 
photon energy of 50 MeV. Figure 1 (left) shows the E tot distribution for the 
data and signal MC near the <fi resonance. One can see good agreement between 
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Fig. 1. The E tot /E cm (left) and x 2 (right) distributions. The points with error bars 
represent experimental events, the histograms show the MC simulation. The arrows 
indicate the cuts imposed. 

the data and signal MC. About 52 x 10 3 events were selected in the whole 
energy range after these requirements. 

Then a kinematic fit requiring energy-momentum conservation (a standard 
4C fit) was performed. We require \ 2 < 15 that provides a good signal/noise 
ratio while the number of rejected signal events is still small (see Fig. 1). 

The reconstruction procedure assumes three photons, i.e., for events with four 
photons a combination of three photons with the minimum \ 2 is chosen. After 
this stage about 48 x 10 3 events remain. 

The dominant background comes from the QED three-photon annihilation: 
e + e~ — > 37. These events can not be completely rejected by selection criteria. 
The 7/7, 7r°7 and background events can be separated using decay dynamics. 
To this end two methods were considered: Dalitz plot analysis and a fit of the 
two-photon invariant mass distributions. 

Figure 2 shows the Dalitz plot for the 37 final state in the meson energy range 
(1011.4 MeV < Eqm < 1027.4 MeV). Here the photons are sorted by their 
energy so that the first photon has a maximum energy: E\ > E^ > E%. The 
Dalitz plot is divided into three regions: D m (340 MeV < E 2 < 385 MeV or 
E 1 < 385 MeV), (491 MeV < E 1 < 511 MeV) and D hg (all the remaining 
events). For each of the three final states (7/7, 7r°7 and QED) we determine 
from the MC simulation the probabilities to enter each region. Based on that, 
from the population of various regions of the Dalitz plot in the data the total 
number of events due to each process is calculated. However, this method can 
provide bias in the signal yield determination because of the possible deviation 
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Fig. 2. The Dalitz plot for the 37 final state at the <t> meson energy. The points 
represent experimental events, the lines indicate boundaries used in the selection 
criteria, see the text for more detail. 

between the signal shape in the data and MC simulation. Additional bias can 
arise from the background processes of the non-QED origin. 

Therefore, we obtain the number of 777 and 7r°7 events by fitting the two- 
photon invariant mass distribution. In this method the signal shape is obtained 
from the data decreasing a possible bias. The difference in the number of se- 
lected events in these two methods (about 3%) was considered as a systematic 
uncertainty because of the separation procedure. 

For 7T° reconstruction the invariant mass of the two softer photons (M23) is 
used. For the 77 signal three combinations are used: 

(1) In case of E\ < ml/y/s, two hard photons are used (M12). 

(2) Otherwise, if E 3 < ml/^/s, we use the first and third photons (M13) . 

(3) In other cases two soft photons (M23) are used. 

Figure 3 shows the two-photon invariant mass distributions for the ir° (left) 
and 77 (right) combinations near the peak of the meson. 

Other possible sources of background are the processes e + e~ — > 777 — > 37r°7, 
e + e~ — ► K^Kl, e + e~ —>■ 77 and e + e~ — ► um° — > 7r°7r°7. The expected number 
of events from these processes was calculated from the detection efficiencies 
determined by the MC simulation and their cross sections independently mea- 
sured at CMD-2 [20,21,22,23,24]. The fraction of background events is negli- 
gible below the <fr meson and is about 2% only in the <fi meson energy range. 
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Fig. 3. The two-photon invariant mass distributions in the 7r° (left) and rj (right) 
mass range. The points with error bars represent experimental events, histograms 
show the MC simulation. 

Above the </> meson the expected cross section of the signal is very low (0.01- 
0.1 nb) and that of the background remaining after all selection criteria has 
close or even higher value. The separation procedure gives 17400 777 events, 
18680 7r°7 events and about 12000 QED events in the whole energy range 
considered. 



3. 1 Approximation of the cross sections 



At each energy point i the cross section of the process o x of a given process is 
calculated using the following formula: 



where N[ is the number of selected events, L x is the integrated luminosity, e- x 
is the detection efficiency and (1 + b\) is the radiative correction at the i-th 
energy point. 

The detection efficiency was calculated from the Monte Carlo simulation tak- 
ing into account corrections obtained from the data and the neutral trigger ef- 
ficiency. The neutral trigger (NT) is part of the CMD-2 trigger system respon- 
sible for events with a final state of photons only, without any charged tracks. 
The NT efficiency was estimated using events of the process e + e _ — > e + e~7 
at each energy point. Its value varied from about 80% to 90%. 
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Fig. 4. The cross section of the process e + e — * r/j. The points with error bars 
represent the experimental data, the curve corresponds to the result of the fit. 

The radiative corrections are calculated according to [25]. The dependence of 
the detection efficiency on the energy of the emitted photon is determined 
from simulation. 

The obtained cross sections of the processes e + e" — > 777, 7r°7 are shown in 
Figs. 4, 5. The observed pattern of the energy dependence is due to the inter- 
ference of the p, u and </> mesons. The detailed information on this analysis is 



7 



Table 1 

The cm. energy, integrated luminosity, number of selected events, detection effi- 
ciency, radiative correction, and Born cross section a of the process e + e~ — ► 777. 



y/s, MeV 


L, nb" 1 


A^cxp 


e, % 


1 + 6 


a, nb 


599.86 


35.2 


0.1 ± 2.3 


2.3 


-0.157 


< 4.71 


629.86 


44.6 


0.7 ±2.8 


7.2 


-0.142 


< 1.68 


659.86 


39.8 


3.0 ±2.9 


11.6 


-0.134 


0.75 ± 0.63 


7f9.86 


56.9 


0.0 ±2.5 


15.0 


-0.126 


< 0.48 


749.86 


42.9 


0.1 ± 1.9 


16.0 


-0.117 


< 0.44 


759.86 


33.7 


3.4 ± 2.7 


15.7 


-0.115 


0.73 ± 0.51 


763.86 


39.7 


5.6 ±3.9 


16.0 


-0.116 


1.00 ±0.61 


769.86 


34.3 


0.2 ±2.1 


17.0 


-0.126 


< 0.63 


773.86 


70.1 


1.0 ± 2.8 


15.8 


-0.147 


0.11 ±0.26 


777.86 


83.6 


2.4 ±3.5 


15.8 


-0.186 


0.22 ±0.27 


779.86 


56.6 


10.0 ±3.9 


14.7 


-0.204 


1.50 ± 0.46 


780.86 


58.5 


2.3 ± 5.9 


8.8 


-0.207 


0.57 ± 1.14 


781.86 


366.8 


41.2 ±9.3 


14.7 


-0.203 


0.94 ±0.18 


782.86 


77.6 


9.8 ±4.4 


15.3 


-0.191 


1.03 ±0.37 


783.86 


71.7 


2.1 ±3.3 


15.7 


-0.172 


0.22 ±0.29 


785.86 


67.0 


8.2 ±3.9 


16.0 


-0.123 


0.87 ±0.36 


789.86 


28.4 


1.5 ±2.1 


15.2 


-0.046 


0.36 ± 0.48 


793.86 


46.2 


1.7 ±2.2 


15.9 


-0.010 


0.23 ±0.31 


799.86 


56.5 


0.9 ±2.0 


15.5 


0.005 


0.10 ±0.23 


809.86 


59.9 


3.9 ± 2.7 


15.9 


0.006 


0.41 ± 0.29 


819.86 


109.4 


8.0 ±4.1 


15.8 


0.004 


0.47 ±0.24 


839.86 


130.4 


8.6 ±4.5 


15.4 


-0.006 


0.43 ± 0.23 


879.86 


167.9 


2.1 ±3.9 


15.2 


-0.037 


0.08 ±0.15 


919.86 


285.4 


5.7 ±5.4 


15.0 


-0.063 


0.14 ±0.13 


939.86 


136.7 


1.1 ±3.8 


15.5 


-0.077 


0.06 ±0.18 


949.86 


226.1 


12.5 ±5.5 


16.2 


-0.085 


0.38 ±0.15 


957.86 


250.1 


6.2 ±4.7 


16.5 


-0.093 


0.17 ±0.12 


969.86 


249.7 


4.7 ±5.1 


17.2 


-0.108 


0.12 ±0.12 


983.93 


307.7 


5.2 ±7.0 


20.4 


-0.132 


0.07 ±0.13 
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Table 2 

The cm. energy, integrated luminosity, number of selected events, detection effi- 
ciency, radiative correction, and Born cross section a of the process e + e~ — ► 777. 



y/s, MeV 


L, nb" 1 


A^cxp 


e, % 


1 + 6 


a, nb 


1003.91 


357.7 


44.0 ± 10.2 


20.3 


-0.192 


0.67 ±0.15 


1010.53 


477.3 


109.6 ± 14.9 


19.9 


-0.227 


1.48 ±0.16 


1015.77 


391.7 


401.3 ±23.0 


20.0 


-0.268 


6.80 ± 0.30 


1016.77 


660.1 


968.1 ±34.8 


19.6 


-0.277 


10.00 ±0.27 


1016.91 


306.1 


497.6 ± 24.9 


20.0 


-0.277 


10.90 ± 0.49 


1017.61 


673.7 


1362.1 ±40.8 


20.1 


-0.282 


13.77 ±0.32 


1017.77 


563.1 


1198.8 ±38.2 


19.9 


-0.282 


14.80 ± 0.44 


1018.58 


410.1 


1230.8 ± 38.4 


20.1 


-0.278 


21.31 ±0.58 


1018.83 


977.5 


2855.9 ± 57.8 


19.9 


-0.274 


21.13 ±0.30 


1019.50 


633.1 


1941.9 ±47.5 


20.1 


-0.254 


21.50 ±0.43 


1019.84 


810.8 


2584.6 ± 54.6 


20.2 


-0.238 


21.54 ±0.73 


1020.62 


876.3 


2231.7 ±51.3 


20.0 


-0.187 


15.52 ±0.30 


1021.54 


440.6 


800.5 ±30.9 


20.0 


-0.112 


9.82 ± 0.36 


1022.79 


551.0 


621.5 ±28.0 


20.1 


0.007 


5.38 ± 0.26 


1027.67 


562.2 


198.9 ± 17.3 


20.2 


0.591 


1.08 ±0.15 


1033.67 


510.8 


100.0 ± 14.0 


19.9 


1.557 


0.38 ±0.14 


1039.59 


447.5 


66.5 ± 11.9 


20.1 


2.911 


0.18 ± 0.13 


1049.80 


312.5 


23.4 ±8.2 


19.6 


6.778 


0.05 ±0.13 


1059.49 


220.6 


9.8 ±5.9 


19.3 


13.272 


< 0.30 


1079.00 


437.0 


4.6 ±6.8 


22.7 


39.899 


< 0.11 


1163.40 


918.2 


0.0 ± 10.3 


21.8 


0.035 


< 0.08 


1310.00 


4249.0 


-0.4 ±21.6 


20.9 


-0.074 


< 0.05 



listed in Tables 1-4. The cross section shown there is a so called "dressed" one, 
which is used in the approximation of the energy dependence with resonances. 
For applications to various dispersion integrals like that for the leading order 
hadronic contribution to the muon anomalous magnetic moment, the "bare" 
cross section should be used [26]. 

The maximum likelihood method is applied to fit the energy dependence of 
the experimental cross sections obtained from the relation (1). 
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Table 3 

The cm. energy, integrated luminosity, number of selected events, detection effi- 
ciency, radiative correction, and Born cross section a of the process e + e~ — > 7r°7. 



v/s, MeV 


L, nb" 1 


N 

1 v cxp 


e, % 


1 + 6 


cr, nb 


599.86 


35.2 


4.8 ±3.6 


12.0 


-0.089 


1.23 ±0.86 


629.86 


44.6 


9.2 ±4.2 


12.8 


-0.093 


1.78 ± 0.74 


659.86 


39.8 


8.6 ±4.4 


12.5 


-0.099 


1.92 ±0.89 


719.86 


56.9 


14.1 ± 5.2 


14.0 


-0.112 


2.00 ±0.65 


749.86 


42.9 


27.1 ±5.9 


14.3 


-0.131 


5.08 ± 0.97 


759.86 


33.7 


35.0 ±6.8 


14.7 


-0.150 


8.31 ± 1.40 


763.86 


39.7 


62.8 ±8.6 


14.5 


-0.162 


12.97 ± 1.52 


769.86 


34.3 


76.7 ±9.2 


15.5 


-0.185 


17.64 ± 1.77 


773.86 


70.1 


281.1 ± 17.3 


14.6 


-0.204 


34.33 ± 1.82 


777.86 


83.6 


721.5 ±27.5 


14.4 


-0.224 


76.65 ± 2.64 


779.86 


56.6 


757.7 ±27.8 


13.8 


-0.229 


125.81 ±4.92 


780.86 


58.5 


717.7 ±27.0 


8.6 


-0.228 


184.98 ± 10.99 


781.86 


366.8 


6619.7 ±82.0 


13.6 


-0.221 


172.26 ±2.47 


782.86 


77.6 


1664.6 ±41.1 


14.9 


-0.206 


183.37 ±4.80 


783.86 


71.7 


1403.6 ± 37.7 


14.9 


-0.183 


162.00 ± 4.72 


785.86 


67.0 


978.8 ±31.6 


13.9 


-0.116 


118.44 ±4.16 


789.86 


28.4 


187.8 ± 13.9 


14.6 


0.050 


42.80 ± 3.80 


793.86 


46.2 


166.2 ± 13.3 


14.8 


0.217 


19.93 ±2.06 


799.86 


56.5 


134.8 ± 12.2 


14.8 


0.441 


11.18 ± 1.52 


809.86 


59.9 


83.7 ±9.9 


15.5 


0.724 


5.22 ± 1.09 


819.86 


109.4 


87.9 ± 10.4 


15.4 


0.906 


2.74 ±0.62 


839.86 


130.4 


61.7 ±9.2 


15.8 


0.901 


1.58 ±0.45 


879.86 


167.9 


17.2 ±6.0 


17.2 


0.342 


0.44 ±0.21 


919.86 


285.4 


20.8 ±6.6 


17.4 


0.021 


0.41 ±0.13 


939.86 


136.7 


18.0 ±5.5 


17.9 


0.001 


0.74 ± 0.22 


949.86 


226.1 


20.1 ±6.2 


18.0 


-0.008 


0.50 ±0.15 


957.86 


250.1 


15.7 ±5.8 


18.4 


-0.015 


0.35 ±0.13 


969.86 


249.7 


11.8 ±5.4 


18.6 


-0.029 


0.26 ±0.12 


983.93 


307.7 


9.4 ±6.3 


19.9 


-0.053 


0.16 ±0.11 
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Table 4 

The cm. energy, integrated luminosity, number of selected events, detection effi- 
ciency, radiative correction, and Born cross section a of the process e + e" 
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4.04 ± U.OZ 


1 090 (59 
1UZU.0Z 


87(5 Q 
0(0.0 


A EL/I 9 _|_ 9(5 O 
404. o It ZD.U 


9n 7 

ZU. ( 


n iQQ 
-u. ioy 


9 8^4-01^ 
Z.00 It U.IO 


1091 KA 
1UZ1.04 


/I /ID (5 
44U.D 


1/19/14-1^9 
140.4 ± 10. Z 


9n fi 

ZU.D 


o n9/i 

-U.UZ4 


1 ^94-017 
1.0Z ± U.I / 


1 099 7Q 
1UZZ. / y 


^1 n 
001. U 


11/1 (5 4-1/1 n 
114. D ± 14. U 


9n a 


n 1 qq 
u.iyy 


n 78 4- n 1 Q. 
U. ( O ± U.lo 


1 097 (57 
1UZ ( .0 1 


^(59 9 
0DZ.Z 


oc 1 4_ Q o 

oo. i ± y.z 


9n fi 

ZU.D 


9 (5(50 
Z.DDU 


(i ns 4- n na 
u.uo ± u.uo 


1 DQQ (57 
lUoo.O ( 


^1 o c 
OlU.o 


1 Q 8 4- 7 £ 

lo.o ± /.0 


9fl 7 

ZU. ( 


A 9 QI (5 
4o.olD 


✓ n 11 
< u.ll 


1 OQO £Q 

luoy.oy 


AA1 K 
44 ( .0 


1(1 /I X fi c 
1U.4 ± 0.0 


9fl 7 

ZU. ( 


79 Q(59 

( z.yoo 


✓ n 11 
< U.ll 


1049 80 


312.5 


1.8 ± 5.1 


20.4 


6.939 


< 0.13 


1059.49 


220.6 


2.0 ±5.3 


20.4 


3.314 


< 0.20 


1079.00 


437.0 


-0.6 ±5.5 


27.4 


1.634 


< 0.08 


1163.40 


918.2 


19.8 ±9.7 


28.1 


-0.047 


0.08 ± 0.04 


1310.00 


4249.0 


48.1 ± 16.2 


26.9 


-0.143 


0.05 ± 0.02 



The Born cross section of these processes can be written as: 

F P7 ( S ) 



(Tp 7 (s) 



V 



A 



V 



\ ay F{m%)' - s - iViTv(s) ' 



(o) 



s 3/2 

m v r v e^ v 



(2) 



where my is the mass of the resonance, r v (s) and T v = r v (m v ) are its 
width at the squared cm. energy s and at the resonance peak (s = m v ), 
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respectively, 5y is its relative phase, F(s) is a factor taking into account the 
energy dependence of the phase space of the final state, i*p 7 (s) = P 7 = — 

mp/2s)) 3 , Oy^ is the cross section at the resonance peak: 



J") 



127T.B 



m 



V^e+e- By^ m 



VJ 



m 



(3) 



v 
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Table 5 



Results of the fits for the 


processes e + e — > 777/ and 


e + e - 




Parameter 


VI 




7T°7 


(°) u 
ay, nb 

<Tu \ nb 

of. ■* 

(°) u 
oj, nb 


0.145 ± 0.063 

n 9qq+ - 175 
u.zyy_ 124 

99 7Q1 +0-220 
la. 1 »i_ .238 

- 




0.70818;^ 

154.821^9 
5.30 ±0.16 

U.±0»_q q 51 


m^, MeV 


1019.52 ±0.05 




1019.46 (fixed) 


rriu, MeV 


782.59 (fixed) 




783.20 ±0.13 


<P</>-wi° 


180 (fixed) 




164.4 ± 7.9 


X 2 /n.d.f. 


72.3/82 




74.0/80 


where 5y^ e +e- and By. 


^ rn are the corresponding 


5 branching ratios. In Eq. (2) 



we sum over all vector mesons relevant at this energy, V = p,u, 0, p', u'. 

The Gounaris-Sakurai model has been used for the description of the p me- 
son [27]. To describe the energy dependence of the uj and meson widths, 
their main decay modes 7r + 7r~7r°, 7r°7 as well as K®K$, K + K~ , 7r + 7T~7r° and 
777, respectively, were taken into account using the same parameterization 
as in [28]. For the p'(1450) the energy dependence of the width assumed 
60% and 40% branching ratios for its decays into ai(1260)7r and ujtt, respec- 
tively [29]. Its mass and width were taken to be 1465 MeV and 400 MeV, 
respectively [12]. The energy dependence of the u/(1420) width is calculated 
assuming the uf — > pn decay. Its mass and width were fixed at the world 
average values of 1425 MeV and 215 MeV, respectively [12]. 

3.2 Results of the fits 

For the fit to the e + e~ — > 777 cross section the resonance cross sections at 
the peak crffl, as well as the meson mass are free parameters. 
The p and u meson phases are chosen to be ip p = (p u = 0° while that for 
the meson is y?^ = 180° in agreement with the quark model. The values of 
the other parameters are taken from Ref. [12]. We also consider a model in 
which in addition to the parameters described above there is a contribution 
tentatively referred to as that of the p'(1450) meson. A fit in the model with 
the p'(1450) doesn't improve x 2 an d results in the value of affi consistent with 
zero, affi = 0.0011q!ooi nl ° an d compatible with our result in the 37r° mode [20] 
af^ = 0.066 ± 0.015 nb. Therefore, for our final results for the 7/7 decay we 
choose a model where a Q , = nb, see Table 5. 
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For the e + e~ — > 7r°7 case the fit parameters are: the cross sections at the reso- 
nance peak uf \ crffi, and the oo meson mass m u . The p — uo phase is fixed 
to the value 13.3° obtained in our study of the process e + e~ — > it + it~ [26]. The 
4> — uj phase is a fit parameter. A fit, which includes a possible ut contribution, 
gives the best \ 2 at the value of a$ significantly differing from zero. Results 
of the best fit are shown in the last column of Table 5. 

3.3 Systematic errors 

There are two types of systematic uncertainties on the cross section <7y: experi- 
mental and model uncertainties. The main sources of experimental systematic 
errors are listed below. The systematic error due to selection criteria is 4% 
estimated by varying the photon energy threshold, total energy deposition, 
total momentum, and x 2 - A possible uncertainty because of the method of 
process separation was estimated to be 3% by comparing our results obtained 
from fitting the distributions of the two-photon invariant mass to those from 
Dalitz plot analysis. The latter method also allows to determine the cross sec- 
tion of the QED process e + e~ — > 37 and it appears to be consistent with the 
theoretical prediction [30]: cr(37) ex p/c"(37)th — 0.973 ±0.018. The uncertainty 
in the determination of the integrated luminosity is 1% and comes from the 
selection criteria of Bhabha events, radiative corrections and calibrations of 
DC and BC. The error of the NT efficiency was estimated to be 2% by trying 
various fitting functions for energy dependence and variations of the cluster 
threshold. The 1% uncertainty of the radiative corrections comes from the 
dependence on the emitted photon energy and the accuracy of the theoretical 
formulae. In total, the experimental systematic uncertainty of the cross section 
is 6%. 

The model uncertainty estimated by comparing the values of the cross section 
at the resonance peak in various models differing by the values of phases 
and resonance parameters was 1% (2%) for the p, 3% (0.1%) for the 00 and 
0.1%(5%) for the meson in the 777 and 7r°7 decay modes, respectively. 



4 Discussion 

In Table 6 we present our results in terms of the product of the branching 
ratios B(V — > e + e~) x B(V — > P7), where P = 77(77°), which is calculated 
from o r y' ) according to (3). For the 777 mode one should additionally take into 
account the branching ratio of the 77 — > 77 decay taking its value B{r\ — > 
77) = (39.43 ± 0.26)% from Ref. [12]. For the 7T°7 mode the corresponding 
value B(ir° -> 77) = (98.798 ±0.032)% from Ref. [12] was included at the MC 
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Table 6 

BjV -> e+e~) x g(y -> P 7 ) 



Decay 


This work 


PDG-2004 


p r/ 7 , 10- 8 


1.50 ±0.65 ±0.09 


1.38 ±0.17 


u; -> 777, 1(T 8 


3.17±i;|f ± 0.21 


3.53 ± 0.35 


<P - r? 7 , 10- 6 


4.093±g;g|§ ± 0.247 


3.85 ± 0.07 


/9 7T° 7 , 10~ 8 


2.90^55 ±0.18 


2.8 ±0.6 


W -► 7T°7, 10~ 6 


6.47 ±0.14 ±0.39 


6 S7+ - 17 
D - J ' -0.15 


7T°7, 10~ 7 


3.75 ±0.11 ±0.29 


3.67 ±0.28 


Table 7 
£(F -► P 7 ) 


Decay 


This work 


PDG-2004 


p -> 777, 10~ 4 


3.21 ±1.39 ±0.20 


3.0 ±0.4 


a; — > 777, 10~ 4 


4.44±f;|| ± 0.28 


4.9 ±0.5 


- r/7, 10- 2 


1.373 ±0.014 ±0.085 


1.295 ±0.025 


p -► 7T°7, 10~ 4 


6.21±^| ± 0.39 


6.0 ± 1.3 


W -► 7T°7, 10~ 2 


9.06 ±0.20 ±0.57 


no+0.28 
°- yz -0.24 


-► 7T° 7 , 10~ 3 


1.258 ±0.037 ±0.077 


1.23 ±0.10 



generation stage. Our results are in good agreement with the world average 
values [12]. 



By dividing the product of the branching ratios above by the corresponding 
world average leptonic width from Ref. [12] one can obtain the branching ratios 
of the radiative decays confronted in Table 7 to the world average values [12]. 

Taking into account a variation of the m u and in various models as well as 
a systematic error caused by the uncertainties of the beam energy calibration, 
we obtain for the resonance masses: 

m u = 783.20 ± 0.13 ± 0.16 MeV, (4) 
m$ = 1019.52 ± 0.05 ± 0.05 MeV, (5) 

consistent with the world average values 782.59 ±0.11 MeV and 1019.456 ± 
0.020 MeV, respectively [12]. 

Our result for the cross section of the process e + e _ — > 777 at the peak of the 
meson can be combined with the independent measurement of the same 
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quantity in the decay mode 77 — > 3n° performed at CMD-2 [20] to obtain the 
ratio of the branching fractions of the 77 meson, B{r\ — > 3tt )/B(i] — > 77). 
Since in both cases the 77 meson decays into neutral particles only, most of 
systematic uncertainties will cancel in such a ratio. As a result, two sources 
of the systematic error survive: 2.5% for the selection criteria and 3% due to 
process separation and we obtain 

K( n _> 377-0) 

f = 0.817 ±0.012 ±0.032, (6) 



which is consistent with the world average value 0.825 ± 0.007 [12]. 



5 Conclusions 

• Using a data sample corresponding to integrated luminosity of 21 pb^" 1 , the 
cross sections of the processes e + e~ — > 777, 7r°7 have been measured in the 
cm. energy range 600-1380 MeV. The following branching ratios have been 
determined: 





-> r? 7 ) = 


(3.21 ±1.39 ±0.20) • IO - 


-4 

j 


B(uj - 


-> 777) = 


(4.44^;| ± 0.28) • 10" 4 , 




m- 


-> 777) = 


(1.373 ±0.014 ±0.085) • 


10 




► 7T°7) = 


(6.21+^1 ±0.39) • 10" 4 , 




B(u - 


► A) = 


(9.06 ±0.20 ±0.57) • 10" 


-2 


B{4>- 


> 7T°7) = 


(1.258 ±0.037 ±0.077) • 


10 



• From the two independent measurements of the — > 777 decay the following 
ratio of the branching ratios of the 77 meson has been obtained: 

B(r] -> 3tt )//3(77 -> 77) = 0.817 ± 0.012 ± 0.032. 
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